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Abstract 

A systematic investigation is given of finite size effects in d = 2 quan- 
tum gravity or equivalently the theory of dynamically triangulated random 
surfaces. For Ising models coupled to random surfaces, finite size effects are 
studied on the one hand by numerical generation of the partition function 
to arbitrary accuracy by a deterministic calculus, and on the other hand by 
an analytic theory based on the singularity analysis of the explicit paramet- 
ric form of the free energy of the corresponding matrix model. Both these 
reveal that the form of the finite size corrections, not surprisingly, depend 
on the string susceptibility. For the general case where the parametric form 
of the matrix model free energy is not explicitly known, it is shown how to 
perform the singularity analysis. All these considerations also apply to other 
observables like susceptibility etc. In the case of the Ising model it is shown 
that the standard Fisher-scaling laws are reproduced. 



A study of finite size effects in statistical systems is of importance for a 
variety of reasons. It is indispensable for a more reliable interpretation of 
numerical data. As shown by Cardy in the context of systems with conformal 
invariance, finite size effects also codify the spectrum of the theory. In most 
cases finite size effects are estimated or parametrised by extrapolating the 
results of simulations carried out for various sizes. It is the purpose of this 
article to show how such finite size effects for two dimensdional quantum 
gravity systems can be studied systematically. 

Let us begin with the case of so called "pure gravity" which is mapped 
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to the one-matrix model described by the partition function 

Z = e- F = Jvcj) exp - Tr[^(f) 2 + V(<j>)] (1) 

where is a iV x iV Hermitean matrix. When V((f>) = g<p 4 ,the expression for 
F in the large-N limit is 

(-12g)"(2n-l)! 
V n\(n + 2)\ 1 ' 

interpreting the coefficient of g n in this sum as the fixed area partition sum 
for an ensemble of random surfaces of fixed area n,one finds 

_ (-12g)"(2n-l)! 
Zn " n!(n + 2)! [6) 

This clearly holds for all n. However, we are interested in the thermodynamic 
or the large-n limit. To this end we use 

e-*3*- 1/2 r 1 1 
FW ^— + Wz + 288F + -I < 4) 

for the Gamma-functionr(2;).This gives 



V47T 8n 

The exponential growth in n is non-universal and depends on V(<f>). The 
power law correction n~ 7 / 2 is universal and defines the string susceptiblity 7 
to be 7 = —6 + 3 when the power law correction is n~ 6 .In this example the 
finite size corrections are of the form [1 — 25 + ..]. 

Next we consider the case of the Ising model coupled to random surfaces. 
This is described by the two matrix model defined by 



Z = e 



= Jv^Vfo exp -Tr^l+^-cfafo + iVM+Vifo))} (6) 



This model is also exactly solvable. But the free energy is not known as an 
explicit function of (c, g). However, a parametric form of the free energy is 
available in the form 



1,0 1 f z dt . . 1 f z dt . N9 1, 9N 3 

i og y _ _ / g ( t ) + / g{ tf + log{l - c 2 ) + - 

2 z g Jo t lg A Jo t 2 4 
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where 

9(z) = (1 * Zz y - c 2 z + 3cV = g (8) 

The numerical generation of Z n is done through the following sequence of 
steps:i)eliminate the logarithms in F by considering instead ii)eliminate 
inverse powers of 1/(1 — 3z) by repeated use of eqn(). These steps result in 
a polynomial of degree 8 for ^ with coefficients depending on (c,g). Again 

repeated use of eqn() written as a quintic in z results in ^ being expressed as 
another quintic. Now the numerical technique to solve this quintic consists 
of writing 

z = 5>i(n)j/ B ; z 2 = J2*2(n)g n ; ;* 5 = £a 5 (n)<7 n (9) 

Then the quintic(eqn()) is equivalent to the recursion relation 

OiH = Cia 1 (n-l)+^2a2(n)+^3a 2 {n-l)+^ 4 a 3 (n)+^ 5 a 4 (n)+^ 6 a 5 (n) (10) 

where £j are coefficients depending on c.It should be noted that all a — i(n) 
vanish for negative n and in fact as z ~ g for small g is the branch we are 
interested in, it follows that a«(n) = for n < ia and that ai(i) = ai(l)\ It 
is also easy to see that ai(l) = (1 — c 2 ) -1 . 

A naive method of solving these recursion relations would make use of 
the identities 

n-1 

a 2 (n) = ai{m) • ai{n — m) 
l 

n-1 

dz{ n ) — J! a l( m ) ' a 2( n — m ) 
1 

n-1 

a 5( n ) — X! a l( m ) ■ a i( n — m ) (11) 
1 

The computational requirements then grow as ~ n 6 if all ai(m);m < n 
are to be evaluated. This is prohibitive. On the other hand the structure of 
the eqns(ll) shows that computation of a 2 (n) only requires the knowledge 
of ai(m) for m < n — 1, computation of a 3 (n) requires only the knowledge 
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ofa2(m) for m < n — 1 and hence of a\{m) for m < n — 2 etc. Thus the eqns 
(10) and (11) can be solved iteratively simultaneously and the comutational 
requirement only grows as ~ n 2 which is entirely manageable. 
The Cosmological Constant Problem 

The main difficulty with the abovementioned method is that aj(n) grow 
very rapidly with n reflecting the (non-universal)exponential growth of Z n .Yoi 
example, Z w ~ 10 100 .In fact by the time one gets to around 100, various 
coefficients exceed the typical machine limits. Of course, special purpose 
arithmetic can be used, but as we shall see soon, a more elegant option is 
available. The idea is to go upto some n m ax such that the coefficients are 
approaching the machine limits and get an estimate for the "cosmological 
constant" jl where Z n ~ e^ n .Then one scales all cij(n) according to 

5i(7i)=e-^ fli(n) (12) 

for all i. This results in the scaling 

6 = e-P £i; £ 3 = e-P £s (13) 

with all other £'s remaining the same. In the statistical mechanical language 
this is equivalent to tuning the system to " criticality" . I will not discuss the 
results here. They can be seen in [|l].That reference also carries the details 
of how a similar analysis can be carried out for the scaling behaviour of 
magnetic susceptibility also. 

Singularity analysis of finite size effects 

This is based on the observation that in a sense n and logg are "conjugate" 
to each other and that the large-n behaviour of Z n controls the convergence of 
the perturbative expansion in g.In particular this is related to the behaviour 
of F near its singular points i.e points where some derivative of F blows up. 

Let us illustrate this first with the example of pure gravity.Here the ex- 
plicit parametric form of F is 

F = -hnz + ^(z-l)(9-z) (14) 



with 



»M = T2?- = » (15) 
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The singular point g c corresponds to zq where g'(zo) = 0. Hence 

z = 2 g c = -1/48 (16) 
Expanding g(z) around Zq one finds 

g-g c - a(z~z ) 2 + b(z-z ) 3 + ... (17) 
Inverting this one gets 

z-z Q ~ c(g-g c ) 1/2 + d(g-g c ) + .. (18) 

Precise values of a,b,c,d can be found in [|IJ. It should be noted that F'(z ) = 
0. Expanding F around z$ one has an expression of the type 

F(z)-F(z ) = A(z-z ) 2 + B(z-z ) 3 + ... (19) 

Combining the last eqns one gets for the singular part of the free energy F 
near g c the expansion 

. . . . 12283^, x ,/ 2 1769472^, x7/2 , v 
F(g) sing ~ F(g c ) = ^-(jg - g c f' 2 + 7— Ms - g c ) 7 ' 2 + .. (20) 

It should be remarked that on the basis of eqns(18,19) one would have ex- 
pected this expansion to start with (g — g c ) 3 ^ 2 .lt turns out that the coefficient 
of this term vanishes. At this stage it appears accidental, but later when we 
present our improved singularity analysis, we will see that such cancellations 
(even more miraculous ones happen in the two-matrix model) have a natural 
explanation. Now we can use the expansion 

(9 ~ 9cT = T.9 n (-9c) a - n r( -Ts (21) 

and the asymptotic expansion of Gamma functions(eqn()) to get the asymp- 
totic expansion for Q n defined by 

(9-9cT = Y.9 n Qn (22) 

n 

giving 

Gnn^oo = (-l) n+1 (-Scr n r(a + 1) n -(l+«)™ e «(l+«)/^ (23 ) 

7T 
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Combining eqns(19,21,22) yields eqn(5). 

In the case of the two-matrix model, the singularity analysis is a bit more 
involved. In the case of the 4 -potential,(? / (zo) = yields z = —1/3 in the 
low temperature phase and a known function Zo(c) in the high temperature 
phase. From the parametric representation it is easy to see that F'(zq) = 
always. Boulatov and Kazakov had claimed that F"(z ) is also zero. But 
we find that this is so only at the Ising critical point c = c cr .Now we present 
the singularity analysis for the cases c ^ c cr and c = c cr separately. 

Non-critical Case 
Here the situation is identical to the pure gravity case except that all the 
coefficients in the expansion depend on c. Again, the coefficient of (g — g c ) 3 ^ 2 
in F, which is now a function of c, vanishes for all clApart from reproducing 
the correct string susceptibility, the form of the finite size corrections one 
obtains are: 

8 76 

HighTempPhase(c = .36) : 1 h ... 

n 

72 69 

LowTempPhase(c = .20) : 1 — + ... (24) 

Critical Case 

In this case we have F"(z ) = 0,g"(zo) = in addition to the vanishing of 
the corresponding first derivatives. Thus eqns(17-19) are now replaced by 

g-g c ~ a(z- z f + b(z- z ) 4 + ... (25) 

z-z ~ c (g-g c ) l l* + d(g-g c ) 2 '' i + h(g-g c ) + .. (26) 

F(z)-F(z ) =A(z-z ) 3 + B(z-Zoy + ... (27) 

Again the (g — g c ) 4 / 3 , (g — g c ) 5//3 terms in the singular part of the free energy 
drop out and the leading singular part ~ (g — g c ) 7 ^ 3 characterstic of the string 
exponent7 = —1/3. The form of the finite size corrections are 

a4 i 0.4287 3.01 1.298 
V(4>) = 4 : 1 + ^7^ iTT- 

„,,, ,o 0.2860 3.05 0.934 

V(4>) = 3 : 1 + ^ 473"- 28 

These eqns represent the most important results of our analysis, namely,that 
the form of the finite size corrections is system-dependent and in this partic- 
ular example depend on the string susceptibility 7. It should be noted that 
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in eqn(28) there are no terms of the type n -2 / 3 . This will become clearer 
with our improved singularity analysis. It should also be noted from eqn(28) 
that the coefficients in the finite size factor are non-universal. This means 
it is possible to choose an appropriate potential that will minimise the finite 
size corrections. This situation is well known in the study of lattice gauge 
theories(" improved action" -principle) . 

On the basis of eqn(28) it had been conjectured in [l[ that if the string 
susceptibility is of the form 7 = p/q with p,q relatively prime, the form of the 
finite size corrections will be 1 H — 57^ + .... We will take up this issue now. 

Improved Singularity Analysis 

The singularity analysis presented above relied on the availability of the 
explicit form of the parametrised free energy. In most cases such explicit 
forms are not available and the above-presented method will fail in those 
cases. There are many models of interest like multi-critical matrix models, 
models with 7 > [H whose finite size analysis one may wish to perform. 

In all such cases where the method of orthogonal polynomials @| can be 
applied, the free energy is erxpressible as 

F = C d£{l-£)lnf{£) (29) 
Jo 

where the function /(£) is defined through the method of orthogonal poly- 
nomials 

J dx w(x)Pi(x)Pj(x) = hidij OneMatrixModels 
J dx dy w(x,y)P i (x)Pj(y) = h^ij TwoMatrix Models (30) 

fi = hi/hi-x ^tv-kx, /(£) (31) 

In what follows we shall only consider one matrix models but consider arbi- 
trary potentials. The analog of the defining eqns for g(z) is now 

g£ = ™(/(0) (32) 

In all these eqns the explicit form of w{x) is dictated by the potential. The 
multicritical points are determined by the vanishing of various derivatives of 
w(x). The upshot of all this is that multicritical models are characterised by 

g£-g 6 = A(f(0 - By + C(/(fl - D)' p+l + .. (33) 
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The meaning of this equation is that all parameters but g have been fixed at 
their critical values and g is near g c . Further, if w(x) is chosen to be the small- 
est polynomial satisfying multicriticality, there is only one term in eqn(38). 

The parametric representation given earlier in terms of the variable z can 
be recovered on identifying /(l) with z whence eqn(37) reads as g = w(z) = 
g(z). Inverting eqn(38) 

f(0=B' +A'(gZ-g c )^ + .. (34) 

On using this with eqn(29) one finds 

m=p— 1 

F sing ~ £ (9-9c) 2+m/p (35) 

m=l 

A naive expectation for the exponent would be 1 + m/p, but the (1 — £) 
factor in the integrand suppresses this dominant exponent. This is basically 
what was behind the surprising cancellations we alluded to earlier. 

Thus we see that in the fractional powers of 1/n characterising finite size 
corrections, two terms will be missing,one corresponding to the first term in 
eqn(35) which controls the leading behaviour and one corresponding to the 
term m = p which is regular. Now we quote results for various cases: 



n 7 Finite Size Correction 



2 


-1/2 


1 + a + 02 + 




3 


-1/3 


1 i , oi _i_ oa 
1 -r ^73 I" ^373 


+ ... 


4 


-1/4 


1 1 ai 1 02 
1 ^ n^Ti "T ^274 


+ ... 



Likewise, the results for models with positive 7 are given below: 
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m 7 Finite Size Correction 



2 


1/2 


1 + 


^ + % + . 




3 


1/3 


1 + 


qi 1 _02 


+ ... 


4 


1/4 


1 + 


a-i i <X2 
^274 "I- ^374 


+ ... 



Indeed there are many interesting issues that need to be properly under- 
stood in this context. For example, in the continuum version of Distler and 
Kawai 0, and of David 0, the liouville field is integrated over all its possible 
values. In such a treatment no "finite size" corrections will ever appear. The 
fixed area formulae hold for all areas. It would be interesting to formulate the 
problem of scaling violations in such approaches. Other issues of both prac- 
tical and theoretical interest are the finite size analysis for other observables 
like the Hausdorff dimension, loop length distributions, resistivity of random 
networks etc. Another outstanding issue is to understand what happens in 
the c = 1 case where 7 = and there are logarthmic modifications to the 
fixed area partition sum. 
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